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We study the nature of the finite-temperature chiral transition in QCD with Nf light quarks in 
the fundamental and adjoint representation. Universality and renormalization-group (RG) argu- 
ments show that the possibility of having a continuous transition is related to the existence of a 
stable fixed point (FP) in the RG flow of a 3D Landau-Ginzburg-Wilson <t> 4 theory with the same 
chiral symmetry-breaking pattern. The RG flow of these theories is studied by field-theoretical 
approaches, computing and analyzing high-order perturbative series, up to six loops. According 
to this RG analysis, the transition in QCD can be continuous only for Nf = 2. In this case it 
belongs to the 3D 0(4) universality class. We also find a stable FP corresponding to a 3D univer- 
sality class with symmetry breaking U(2)z,®U(2)s^U(2)y, which implies that the transition can 
be continuous also if the axial-anomaly effects are suppressed at T c . In the case of quarks in the 
adjoint representation, we can have a continuous transition for Nf = 1,2. For Nf = 1 it belongs 
to the 0(3) universality class. For Nf = 2 it belongs to a new 3D universality class characterized 
by the symmetry breaking SU(4)^SO(4). 
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The thermodynamics of Quantum Chromodynamics (QCD) is characterized by a transition 
from a low-J hadronic phase, in which chiral symmetry is broken, to a high-J phase with de- 
confined quarks and gluons (quark-gluon plasma), in which chiral symmetry is restored [j|]. Our 
understanding of the fmite-r phase transition is essentially based on the relevant symmetry and 
symmetry-breaking pattern (SBP). In the presence of Nf light quarks the relevant symmetry is the 
chiral symmetry U(l)v <8> SU(iV/)i SU(Nf)g. At T = this symmetry is spontaneously broken 
to U(l)y<8>SU(Af/)v with a nonzero quark condensate The finite-r transition is related to the 

restoring of the chiral symmetry. It is therefore characterized by an Nf x Nf complex-matrix order 
parameter <&,,-, related to the bilinear quark operator YuYRj> and the SBP 

SU(Nf)L®SU(Nf)R^>SU(Nf) v . (1) 
If the axial U(1)a symmetry is effectively restored at T c , the expected SBP becomes 

U(N f ) L ®\J(Nf) R -» U(N f ) v . (2) 

Lattice Monte Carlo (MC) simulations suggest that these is not the case in three-color QCD. How- 
ever, since the anomaly gets suppressed in the large-A^ c limit (d^J^ ^-Q), SBP (||) may be relevant 
in the large-A^ limit. 

Although deconfinement and chiral symmetry restoration are apparently related to different 
mechanisms, they seem to be somehow coupled in QCD, since lattice computations show that 
the Polyakov loop has a sharp increase at T c where the chiral condensate vanishes. However, the 
interplay between the two effects is not clear yet. Insight into this question may be gained by 
investigating related models, such as S\J(N C ■) gauge theories with Nf Dirac fermions in the adjoint 
representation (aQCD). 1 Unlike QCD, aQCD is also invariant under global TL^ c transformations 
related to the center of the gauge group S\J(N C ), as in pure SU(iV c ) gauge theories. There are 
two well-defined order parameters in the light-quark regime, related to the confining and chiral 
modes, i.e. the Polyakov loop and the quark condensate. One generally expects two transitions: 
a deconfinement transition at associated with the breaking of the TL^ c symmetry, and a chiral 
transition at T c in which chiral symmetry is restored. In aQCD with Nf massless flavors the chiral- 
symmetry group extends to SU(2Afy). At T = this symmetry is expected to spontaneously break 
to SO(2Nf), due to quark condensation. Therefore the SBP at the finite-r chiral transition is 
expected to be 

SU(2tf/) -» SO(2N f ) (3) 

with a symmetric 2Nf x 2Nf complex matrix as order parameter related to the bilinear quark con- 
densate. If the axial \J(1)a symmetry is restored at T c , the SBP is 

\]{2Nf)^0{2N f ). (4) 

MC simulations for N c = 3 and Nf = 2 |3j] show that the deconfinement transition at Tj is first 
order. Results at the chiral transition appear consistent with a continuous transition. Interestingly, 
the ratio between the two critical temperatures turns out to be quite large, T c /T^ pa 8, suggesting a 
rather weak interplay between the corresponding underlying mechanisms. 

1 AQCD is asymptotically free only for Nf < 1 1 /4, thus only the cases Nf = 1,2 are interesting. 
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In order to study the nature of the finite-J chiral transition in QCD and aQCD, we exploit uni- 
versality and renormalization-group (RG) arguments, as originally applied by Pisarski and Wilczek 
in Ref. jQ). They can be summarized as follows. 

(i) Let us first assume that the phase transition at T c is continuous for vanishing quark masses. In 
this case the length scale of the critical modes diverges approaching T c , becoming eventually much 
larger than Y/T c , which is the size of the euclidean "temporal" dimension at T c . Therefore, the 
asymptotic critical behavior must be associated with a 3D universality class with the same SBR 

(ii) The existence of such a 3D universality class can be investigated by considering the most gen- 
eral Landau-Ginzburg-Wilson (LGW) <I> 4 theory compatible with the given SBR which describes 
the critical modes at T c . Neglecting the U(1)a anomaly, it is given by 



If <J>,y is a generic N xN complex matrix, the symmetry is \J(N)l®U(N)r, which breaks to \J(N)v 
if vo > 0, thus providing the LGW theory relevant for QCD with Nf = N. If <1> !; is also symmetric, 
the global symmetry is \J(N), which breaks to 0(/V) if vq > 0, which is the case relevant for aQCD 
with Nf = N/2. The reduction of the symmetry to SU(iV/) L <g>SU(Ar/) R for QCD [SV(2N f ) for 
aQCD], due to the axial anomaly, is achieved by adding determinant terms, such as 



Nonvanishing quark masses can be accounted for by adding an external-field term Tr (//<!> + h.c). 
(iii) The critical behavior at a continuous transition is determined by the fixed points (FPs) of the 
RG flow: the absence of a stable FP generally implies first-order transitions. Therefore, a necessary 
condition of consistency with the initial hypothesis (i) of a continuous transition is the existence of 
stable FP in the corresponding LGW <I> 4 theory. If no stable FP exists, the finite- T chiral transition 
of QCD (aQCD) is predicted to be first order. If a stable FP exists, the transition can be continuous, 
and its universal critical behavior is determined by the FP; but this does not exclude a first-order 
transition if the system is outside the attraction domain of the stable FP. 

The above arguments show that the nature of the finite-r transition in QCD and aQCD can be 
investigated by studying the RG flow of the corresponding 3D LGW <J> 4 theories. For this purpose 
we consider two different field-theoretical (FT) perturbative approaches. One is defined within the 
massive (disordered) phase using a zero-momentum renormalization (MZM) scheme. The other 
one is defined within the massless (critical) theory using a minimal subtraction (MS) scheme (we 
consider a 3D-MS scheme without e expansion). The RG flow is determined by the FPs, which 
are given by the common zeroes of the /3 -functions associated with the quartic couplings. A FP 
is stable if all eigenvalues of its stability matrix have positive real part. Using symbolic manipula- 
tions programs, we computed the expansion up to six loops in the MZM scheme (which requires 
the calculation of approximately 1500 Feynman diagrams) and up to five loops in the MS scheme. 
Since perturbative FT expansions are asymptotic, it is necessary to resum the series. This is done 
by exploiting Borel summability and knowledge of the large-order behavior, which is inferred by 
semiclassical calculations of instanton solutions. We first resum the /3 -functions, and then search 
for their common zeroes. This computation is essentially nonperturbative, because the resumma- 
tion uses nonperturbative information on their large-order behavior. The comparison of the analyses 



jgf u(JV) =Tr(d^)(d^) + rT^+'^ (TrO^) 2 + ^Tr (O 1 ^) 2 . 



(5) 




(6) 
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of the MZM and 3D-MS expansions provides nontrivial crosschecks of the results. Details of these 
calculations can be found in Refs. ^ ^] (see also [Q]). In the following we summarize the main 
results. 

For N =\, the case relevant to Nf = 1 QCD, J^u(i) reduces to the 0(2) symmetric <1> 4 theory, 
corresponding to the XY universality class of superfluid transition in 4 He, see, e.g., The de- 
terminant term related to the axial anomaly, cf. Eq. (^), plays the role of an external field, thus no 
continuous transition is expected, but a crossover. 

The case N = 2 without anomaly, cf. Eq. (Q), relevant for Nf = 2 QCD, was originally analyzed 
by Pisarski and Wilczek [Q] within the e = 4 — D expansion to one-loop order, which means close 
to 4D. No stable FP is found close to 4D, as also in the case of symmetric matrix field relevant for 
Nf = 1 aQCD. Thus a naive extrapolation to 3D indicates first-order transitions. However, in some 
physically interesting cases the extrapolation of £ -expansion calculations to e = 1 fails to provide 
the correct physical picture in 3D: for example, in the Ginzburg-Landau model of superconductors 
[10], and in 0(2)(g>0(/V) theories of some frustrated spin models [|l~i|]. Actually, this happens also 
in the case of the LGW <J> 4 theory @ with N = 2. Indeed, the analysis of the high-order series 
shows the presence of stable FPs in both MZM and 3D-MS schemes, contradicting earlier analyses 
based on the e expansion around 4D. 2 These results imply the existence of 3D universality classes 
with SBP U(2) L 0U(2) s ^U(2)v and U(2)^0(2), 3 corresponding respectively to Nf = 2 QCD and 
Nf = 1 aQCD in the case of suppressed \J(1)a anomaly. 

In two-flavor QCD, taking into account the U(1)a anomaly, SBP (|l|) becomes equivalent to the 
one of the 0(4) vector universality class, i.e. SO(4)^SO(3). SBP (g) of N f = 1 aQCD is instead 
equivalent to the one of the 0(3) vector universality class. This means that, if the transition is 
continuous, it must show the 0(4) scaling behavior in Nf = 2 QCD and the 0(3) one in Nf = 1 
aQCD. See, e.g., Ref. [§] for a recent review on the 0(/V) universality classes. Actually, the LGW 
<I> 4 theory corresponding to these cases is 

-2su(2) =Tr(^<D t )(^<I>) + rTr^ + ^ (Tr^<I>) 2 + ^Tr (O^) 2 (7) 

+w (detO 1 + detO) + j (TrO 1 *) (detO t + detO) + ^ [(detO 1 ") 2 + (det<I>) 2 ] , 

where Wo,xo,yo ~ g and g parametrizes the effective breaking of the U(1)a symmetry. If the 
anomaly is suppressed (g = 0), then wo = xq = yo = 0. -£?su(2) contains two quadratic (mass) 
terms, therefore it describes several transition lines in the T-g plane, which meet at a multicritical 
point for g = 0. In the case of QCD the multicritical behavior is controlled by the \J(2)l<8>\J(2)r 
symmetric theory. Possible phase diagrams in the T-g plane are shown in Fig. |l[ When g ^ the 
transition may be first order or continuous in the 0(4) universality class. Actually, we may also 
have a mean-field behavior (apart from logarithms) for particular values of g, see Fig. [j] A similar 
scenario applies also to Nf = 1 aQCD. 

No stable FPs are found for Af > 2 in the LGW theory (||). Thus, neglecting the anomaly, 
transitions are always first order when Nf > 2 for QCD and Nf > 1 for aQCD. In most cases this 
result does not change if we take into account the axial anomaly, cf. Eq. (Q). The only exception is 
the case related to the Nf = 2 aQCD, where we find a stable FP corresponding to a 3D SU(4)/SO(4) 

2 This result was overlooked in Ref. Some details can be found in Ref. [Q], 

3 The U(2)/0(2) universality class is also the one of the normal-to-planar superfluid transition in 3 He [fla 
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Figure 1: Possible phase diagrams in the T-g plane for the LGW theory ([7|) describing the transition of 
Nf = 2 QCD, in the case the transition at the multicritical point, i.e. for g = 0, is continuous (left) or first 
order (right). Thick black lines indicate first-order transitions. At their end points, thus for particular values 
of g, the transition should be of mean-field type (apart from logarithms). 





U(1)a anomaly 


suppressed anomaly at T c 


QCD 


SU(N f ) L ® S\J(N f ) R -> S\J{N f ) v 


U(Nf) L ®\J(N f ) R ^\J(Nf) v 


N f =l 


crossover or first order 


0(2) or first order 


N f = 2 


0(4) or first order 


U(2) L ®U(2) /? /U(2)y or first order 


N f > 3 


first order 


first order 


aQCD 


SU(2A^ / ) -» SO{2N f ) 


U(2N f ) -» 0(2N f ) 


N f = l 


0(3) or first order 


U(2)/0(2) or first order 


N f = 2 


SU(4)/SO(4) or first order 


first order 



Table 1: Summary of the RG predictions. For each case we report the possible types of transition, indicating 
the universality class when the transition can also be continuous. 



universality class, with critical exponents v « 1.1 and v\ « 0.2. Note that, although SU(4)~0(6), 
the SBP of Nf = 2 aQCD, i.e. SU(4)^SO(4), differs from that of the 0(6) vector model, i.e. 
SO(6)^ SO(5). Thus, according to the standard paradigm that relates the universality class to the 
SBP, the corresponding universality classes and critical behaviors must be different. 

The predictions of our RG analysis for the finite-r chiral transitions in QCD and aQCD 
are summarized in Table |l[ These transitions have also been investigated by lattice MC simula- 
tions. Overall, MC results for two-flavor QCD seem to favor a continuous transition in the con- 
tinuum limit. However, a satisfactory check of the 0(4) scaling behavior has not been achieved 
yet. Results obtained using Wilson fermions appear consistent with a continuous transition in the 



0(4) universality class [|13|]; MC simulations using staggered fermions appear more problematic 
pi, |16[ H7), |||. Unlike Nf = 2 QCD, the transition scenario appears settled for N f >3: MC 
simulations [TBI |18|] show first-order transitions, in agreement with the RG predictions. Finally, in 
the case of Nf = 2 aQCD the available MC results [Q, ||] favor a continuous transitions. But they 
are not yet sufficiently precise to check the critical behavior of the 3D SU(4)/SO(4) universality 
class. 



199/5 



Finite-temperature chiral transition in QCD with quarks in the fundamental and adjoint representation 
Ettore Vicari 



References 

[1] See, e.g., F. Wilczek, QCD in extreme conditions, hep-ph/0003183; F. Karsch, Lectures on Quark 
Matter, hep-lat/010 6019. 

[2] F. Karsch and M. Lutgemeier, Deconfinement and chiral symmetry restoration in an SU(3) gauge 
theory with adjoint fermions, Nucl. Phys. B550 (1999) 449 [hep-lat/9812023]. 

[3] J. Engels, S. Holtmann, and T. Schulze, Scaling and Goldstone effects in a QCD with two flavours of 
adjoint quarks, Nucl. Phys. B724 (2005) 357 [hep-lat/ 0505008]; The chiral transition ofNf = 2 
QCD with fundamental and adjoint fermions, PoS (LAT2005) 148 [hep-lat / 5 9 1 0]. 

[4] R.D. Pisarski and F. Wilczek, Remarks on the chiral phase transition in chromodynamics, Phys. Rev. 
D29 (1984) 338. 

[5] F. Basile, A. Pelissetto, and E. Vicari, The finite-temperature chiral transition in QCD with adjoint 
fermions, JHEP 02 (2005) 044 [hep- 1 h / 4 1 2 2 ] . 

[6] A. Butti, A. Pelissetto, and E. Vicari, On the nature of the finite-temperature transition in QCD, JHEP 
08 (2003) 029 [hep-ph/0307 03 6]. 

[7] P. Calabrese and P. Parruccini, Five-loop e expansion for U (n) ® U(m) models: finite-temperature 
phase transition in light QCD, JHEP 05 (2004) 018 [hep-ph/ 4 0314 0]. 

[8] M. Campostrini, M. Hasenbusch, A. Pelissetto, P. Rossi, and E. Vicari, Critical behavior of the 
three-dimensional XY universality class, Phys. Rev. B 63 (2001) 214503 [cond-mat / 00 103 60]. 

[9] A. Pelissetto and E. Vicari, Critical phenomena and renormalization-group theory, Phys. Rep. 368 
(2002) 549 [cond-mat/0012164]. 

[10] S. Mo, J. Hove, and A. Sudb0, Order of the metal-to-superconductor transition, Phys. Rev. B65 
(2002) 104501 [cond-mat /010 92 60]. 

[11] P. Calabrese, P. Parruccini, A. Pelissetto, and E. Vicari, Critical behavior of 0(2)®0(N) symmetric 
models, Phys. Rev. B70 (2004) 174439 [cond-mat/0405667]. 

[12] M. De Prato, A. Pelissetto, and E. Vicari, Normal-to-planar superfluid transition in 3 He, Phys. Rev. 
B70 (2004)214519 [cond-mat/0312362]. 

[13] A. Ali Khan et al. (CP-PACS Collaboration), Phase structure and critical temperature of two-flavor 
QCD with renormalization-group improved gauge action and clover improved Wilson action, Phys. 
Rev. D63 (2001)034502 [hep-lat /0008011]. 

[14] F. Karsch, E. Laermann, and A. Peikert, Quark mass and flavor dependence of the QCD phase 
transition, Nucl. Phys. B605 (2001)579 [hep-lat/0012023]. 

[15] J. B. Kogut and D. K. Sinclair, Scaling behavior at the N t — 6 chiral phase transition for 2-flavor 

lattice QCD with massless staggered quarks and an irrelevant 4-fermion interaction, Phys. Rev. D64 
(2001)034508 [hep-lat /0104011]. 

[16] J. Engels, S. Holtmann, T. Mendes, and T. Schulze, Finite-size-scaling behavior for 3d 0(4) and 0(2) 
spin models and QCD, Phys. Lett. B514 (2001) 299 [hep-lat/0105028]. 

[17] M. D'Elia, A. Di Giacomo, and C. Pica, Two flavor QCD and Confinement, hep-lat / 0503030. 

[18] Y. Iwasaki, K. Kanaya, S. Sakai, and T. Yoshie, Chiral phase transition in lattice QCD with Wilson 
quarts, Z. Physik C71 (1996)337 [hep-lat/9504019]. 



199/6 



